Aharonov-Bohm magnetization of mesoscopic rings caused by inelastic relaxation. 
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The magnetization of a system of many mesoscopic rings under non-equilibrium conditions is 
considered. The corresponding disorder- averaged current in a ring I{cf>) is shown to be a sum of the 
'thermodynamic' and 'kinetic' contributions both resulting from the electron-electron interaction. 
The thermodynamic part can be expressed through the diagonal matrix elements J^^ of the current 
operator in the basis of exact many-body eigenstates and is a generalization of the equilibrium 
persistent current. The novel kinetic part is present only out of equilibrium and is governed by 
the off-diagonal matrix elements J^^. It has drastically different temperature and magnetic field 
behavior. 
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Recently there has been a considerable interest in 
mesoscopic phenomena under non-equilibrium condi- 
tions. The relaxation of electron energy distribution 
function due to inelastic processes of electron interac- 
tion has been measured in mesoscopic wires Q]. The 
Kondo effect in quantum dots with finite bias has been 
intensively discussed theoretically j^, ^ and observed ex- 
perimentally [Q. 

The non-equilibrium effects are also important for the 
magnetic response of a system of mesoscopic rings 
^. In equilibrium a weak constant magnetic field gives 
rise 1^] to a magnetization of a system of many normal- 
metal mesoscopic rings that corresponds to the average 
persistent current per ring 



d(j) 



Y^I^r^'^ sin 



(1) 



where is the magnetic flux through a ring, is the 
total energy of interacting electrons in a ring averaged 
over the disorder and the thermal ensembles, and 0o — 
hc/e is the flux quantum. 

The remarkable feature of the flux dependence in 
Eq.(|l]) is its periodicity and odd character. The latter is 
due to the dissipationless nature of the equilibrium per- 
sistent current. Indeed, in equilibrium both directions of 
time are equivalent, so that time-reversal symmetry re- 
quires /(^'-^^((/)) = — /^^'-^'(— 0). The periodic and odd in 
(/) magnetization of the form Eq.(^ has been indeed ob- 
served in a sample containing 10^ mesoscopic copper 
rings. 

However the disorder-averaged current in a ring must 
respect another symmetry related with the space homo- 
geneity of the disorder-averaged system. This is space 
reflection about the ring diameter. It is easy to see that 
both the current and the magnetic flux change sign un- 
der such a reflection, so that the symmetry relationship 



= should hold even in the case where equi- 

librium and time reversal symmetry is not assumed. We 
thus arrive at the statement that the odd character of 
7(0) cannot be used as an evidence that the disorder- 
averaged current observed in Q is an equilibrium persis- 
tent current. The same is true for the 0o/2 periodicity. 

It has been shown by straightforward calculations |6| 
that a non- equilibrium dc current I^'^'^\<j)) of the same 
form as Eq.(|l|) indeed arises when the ring is driven out 
of equilibrium by an external ac electric field. 

A non-equilibrium current is sustainable for a reason- 
ably long time only if there is an external force acting 
on the electron system. This force can produce a non- 
equilibrium dc current either by a direct action on elec- 
trons in the ring or indirectly. The case of the external 
ac electric field in Refs.[|[ |j can be considered as a rep- 
resentative example of the direct effect. 

In the present Letter we consider another, indirect 
mechanism of a non-equilibrium current which is related 
with the relaxation of the given non- equilibrium electron 
energy distribution created by an external force. The 
similar relaxation-induced mesoscopic photovoltaic effect 
has been considered in Ref.||]. However, in contrast to 
Ref. 1^ we do not assume an electron-phonon mechanism 
of relaxation which is extremely weak at low tempera- 
tures [ p^ . Instead we consider much more effective in- 
elastic processes due to electron- electron interaction. 

One can imagine the experimental geometry where the 
direct action of the external force on electrons in the ring 
is absent while the one-particle energy distribution devi- 
ates from the Fermi-Dirac form. Consider a mesoscopic 
ring weakly coupled to the center of a wire connecting two 
reservoirs with the temperature T and the chemical po- 
tential difference V (see Fig. la) maintained by an applied 
voltage. For a sufficiently short wire the distribution of 
one-particle energies in the ring f{E) is roughly a super- 
position of two Fermi-Dirac distributions corresponding 
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to two reservoirs |^ (see Fig. lb). 
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FIG. 1: a). The experimental geometry, b). Non-equilibrium 
energy distribution function with two steps, c). DoS for a 
Kondo system 

Below we identify the disorder-averaged do current 
/^''^ {(f)) in a ring which is entirely due to the relaxation 
of this distribution resulting from the electron- electron 
interaction. We show that this relaxation-induced cur- 
rent can be represented in the form Eq.(|l|) with the flux 

(r) 

harmonics In expressed through the electron energy dis- 
tribution function f{E) = i(l — He)- 



T(r) _ (-1 



hg 



dEdE' 
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dRUE',E) 
doj 



(2) 

where — e is the electron charge, g is the dimensionless 
conductance of the ring with the perimeter L and the 
cross-section area S, C„ = ■(2^^m=i 
and 
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RUE,E') = [{hE^hE-^){l-hE'hE'-^) (3) 
- {hE' - hE'-c^)il ~ hshE-uj)]- 

The diffusion coefficient De = Dq -\- SDe is supposed 
to have a small but essentially energy-dependent part 
SDe <Si Dq. It is also assumed that the characteristic 
energy scale of the energy dependence in 5De is larger 
than the typical energy transfer uj ~ Et — TlDq/L'^. 

The function R^{E,E') in Eq.(^ is exactly the com- 
bination of the electron energy distribution function h e 
that enters in the inelastic collision integral 



St{E)^ J dE'dwP{uj)R^{E,E'). 



(4) 



For the Fermi-Dirac distribution He ^ tanh(_E/2r) we 
have identically Rui{E,E') = 0, and both the relaxation 
rate Eq.(|^) and the non-equilibrium curr ent Eq.(|) van- 
ish. This quantifies an intimate relationship between 
them. 

Eqs.(||),(^) is the main result of the paper. It is valid 
for the case of a pure potential disorder with no spin-orbit 
interaction (orthogonal symmetry class) for g ^ 1 as long 
asT :^ Et and n < L^/L, where = (Det!^™^^)^ > 



L is the dephasing length and rjp^'^^^ is the dephasing 
time in the ring. 

The non- Fermi-Dirac form of He is the necessary but 

(r) 

not the sufficient condition for In to be non-zero. The 
global balance in Eq.(^) requires J dE St{E) = which 
results from the obvious identity / dEdE' Ruj{E, E') — 0. 



Then one concludes from Eq.(H) that /„ 



(r) 



unless 



the electron diffusion coefficient De — v'^te/3 is energy- 
dependent. 

For the white-noise impurity potential within the non- 
crossing (self-consistent Born) approximation the prod- 
uct of the elastic scattering time te and the disorder- 
averaged one-electron density of states (DoS) ve is en- 
ergy independent teve = const. If the energy depen- 
dence of i^E is stronger than the dependence of the elec- 
tron velocity v{E), we obtain De = D^v/ve, where v 
is the (one-spin) DoS outside the region (near the Fermi 
energy) of the strong energy dependence. 

Though the nature of the ^-dependence of ve goes 
beyond the scope of this paper we note that the electron- 
electron interaction results in the strong energy depen- 
dence of the one-particle tunnel density of states j/^""-* 
exactly at the Fermi level Q ^ 0. In our problem 
the corresponding DoS i^^'' should be obtained from the 
non-perturbative treatment of the electron-electron inter- 
action and will be considered elsewhere. An alternative 
mechanism of the energy dependence (see Fig.lc) is the 
Abrikosov-Suhl peaks in ve &t E = ±y/2 that arise 
because of the Kondo effect. 

Now we outline the derivation of Eq.(||). We start 
with the expression for current density J in terms of 
the components of the matrix Green's function G = 
G^' 
G^ 



in the Keldysh technique ||T^, |l^ 



J^-Tr{JG^]^ 



(5) 



where J is the current density operator, G^'^ are the 
retarded (advanced) electron Green's functions, and 
is the Keldysh function. The matrix Green's function 
G = Gq + iGq Go + iGjQ Gg is calculated in 
the first order in the screened electron interaction U_ — 

Q jjA 1 1 where Gg is the matrix Green's function 

without electron interaction in the presence of a static 
disorder potential and T^^'^ are the Fock and the Hartree 
self-energy parts due to electron interaction. 

The electron energy distribution function He enters in 
the theory through the Keldysh component of the unper- 
turbed matrix Green's function Gg via the ansatz [pi] , [l6| 



G^{E)^hE [G^iE)-G^{E)]. 



(6) 



The effective interaction IJ_ is calculated in the random- 
phase approximation which after averaging over disorder 
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in the limit g ^ 1 yields: 

([/fl _C/^)^ q = j dE'U^,^{E'){hE' -hE'-^), (7) 

wheie2ivDoq^U^,^{E') = {Dlq^ +uj^)/{Dl,q^ + lu^). It 
reduces to the well known expression in the case of an 
energy-independent diffusion coefficient SDe' = if one 
uses the identity / dE' {Iie' —fiE'-Lu) = The Keldysh 
component of the effective interaction U;^^ is given by 
Eq.(^ where {He' —hE'-u) is replaced by {l — hE'tiE'-ui)- 
Rewriting the self-energy parts E^'^ explicitly in 
terms of the components G^'^'^ and using the ansatz 
Eq.(^ one identifies three different interaction- 

induced contributions to the current Eq.(^: 

11 = {AAR- RATl)[ihE ~ hE-^)U^ - (8) 
- {1 - hEhE-^) {U^ ~ U^)] 

12 = ARIi{US-pUi){l^hEhE-u) 

+ RAA {U^ - pU^) (1 - hEhE-o.) 

13 = {RRR- AAA)hEU^ - (l-hEhE^u.) X 

[RRR ([/« - pU^) + AAA (U^ - pU^)], 

where AAR = G^{E - uj)G^{E)JGf{E), and the inte- 
gral over all E and w is assumed in Eq.(||). The degener- 
acy factor p = 2 for unpolarized electron spins and p = 1 
if spins are fully polarized by the parallel magnetic field. 

The three contributions are very different in charac- 
ter. Since I3 contains only retarded or only advanced 
Green's functions it vanishes after averaging over dis- 
order at a constant chemical potential. The contribu- 

(PC) 

tion I2 amounts to /A . In particular, the result of 
Ref . follows from this contribution if one assumes the 
equilibrium Fcrmi-Dirac distribution and the multiplic- 
ity p ~ 2. Indeed, using the identity L-^Gq'^{E) = 

-G^'^{E)JG^'^{E) one concludes that I2 can be ex- 
pressed in terms of the flux derivative of the effective en- 
ergy functional £cs{4>) that in the absence of equilibrium 
stands for the total energy f (0) in Eq.(|l|). In order to 
understand the physical meaning of the additional contri- 
bution Ii we invoke the basis of exact many-body electron 
states ^'^ and the corresponding matrix elements of the 
current operator J^^. In this representation the 'ther- 
modynamic' contribution I2 is expressed in terms of the 
diagonal matrix elements J^^ only. In contrast to that 
the 'kinetic' contribution Ii contains only off-diagonal 
matrix elements J^^. One can check that the replacing 
the current operator J in the expression for Ii by the unit 
operator results in vanishing of the whole expression, as 
it is required by the orthogonality of the different many- 
body wave functions. Thus the contribution Ii contains 
an information on the overlap of the different many-body 
wave functions that is totally absent in I2. 

We also note that Ii comes entirely from the Fock-type 
diagrams, while I2 contains both the Hartree and the 



Fock contributions. The two contributions have oppo- 
site signs and the balance between them is controlled by 
the multiplicity factor p. Therefore the thermodynamic 
contribution I2 is very sensitive to the parallel magnetic 
field that leads to the cancellation of the main part of 
the persistent current which is due to the real part of the 
screened interaction U^'^. At the same time the effect of 
the parallel magnetic field on the kinetic contribution re- 
duces merely to a possible variation of the -E-dependence 
of DoS. 

The contribution (Ii) averaged over disorder is exactly 
the relaxation- induced non-equilibrium current /*•''-'((/)). 
Substituting Eq.(^ into the first of Eqs.(^) we obtain 

-— ^ / du;J„{E,E',u;)RUE,E'), 

where R^{E,E') is given by Eq.(|) and J„{E,E',lu) is 
the n-th flux harmonic of 

- J2iiG^ - G^)e^. G^{E)JG^{E)), U^,^{E'). (10) 

The disorder average (...) in Eq.([l0|) is done within the 
impurity diagrammatic technique. In the leading ap- 
proximation in 1/g the result is described by the di- 
agram Fig. 2. The quasi-one-dimensional geometry of 



k + q,0 




FIG. 2: Disorder averaging of the ylAR term in Eq.(|lO|). 
Solid lines labelled by R{A) axe disorder averaged electron 
Green's functions {Gq'^^\E)) with the energy indicated; the 
wavy lines are the diffusion propagators (diffusons and cooper- 
ons) with the momentum and frequency indicated; the dashed 
square is the Hikami box with the vector vertex J denoted by 
the bold triangle; the bold solid line is the screened electron 
interaction Uui.qiE'). 

rings is taken into account by the quantization of mo- 
menta of diffusion propagators [Deci^ ~ iuj]^'^ and elec- 
tron interaction U^_q(E'), such that the transverse mo- 
mentum q± — 0, and the longitudinal momentum q|| = 
(27r/L) (m — 2(/)/(/)o) for the cooperons and g|| = q = 
{2'Km/L) for the diffusons and the electron interaction, 
where m = 0, ±1,±2,... is an integer. The only con- 
straint is that the electron interaction must be zero at 
zero momentum. Performing the Poisson summation 
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over k, we obtain 



(1 



q=iO 



(11) 

where = [De / {—iio)]^ and g — vDS/L. 

The energy dependence of Jn{E,E',uj) originates i). 
from the i?-dependence of the triangle of Green's func- 
tions in Eq.(^o|) and ii). from the iJ'-dependence of the 
polarization bubble in the effective interation Ui^,q{E'). 
It is easy to see that only the latter is important. Indeed, 
let us expand R^{E,E') in Eq.(H) up to the linear in uj 
term. If I?^;' = Dq one can perform the Wick rotation 
—iuj UJ that immediately gives Jn{E, E', to) toduj — 
because of the 3 sign in front of the sum in Eq.(^l]). 
This is not true if SDe' / 0, as in this case the Wick 
rotation leads to divergency. A careful analysis shows 



that the term proportional to SDe in Eq.(ll) makes a 
contribution to Eq.(H) that is small by the parameter 
[Et/TY^'^ ^ 1 compared to that resulting from SDe'- 
Neglecting this contribution by setting De — Dq we ar- 
rive at a finite result Eq.lQ). 

Note that Jn{E,E',uj) is not exponentially small at 
oj ^ Et or L ^ L^. This can be traced back to 
the structure Gq '^{E - uj)Gf {E) J (E) of the kinetic 
term Ii that allows to build a cooperon at zero frequency 
|lo| . This is impossible for the thermodynamic term I2 in 
which both Green's functions with the same energy E are 
either retarded or advanced. Therefore the correspond- 
ing kernel for the thermodynamic term is proportional to 
exp[— L/L;j] ~ exY'l—yjT I Et\- This is the reason why 
the kinetic term wins over the thermodynamic one for 
T:^ Er- 
in conclusion, using the Keldysh formalism we identi- 
fied two different contributions, thermodynamic and ki- 
netic, to the disorder-averaged magnetization of meso- 
scopic rings with a non-equilibrium distribution of one- 
electron energies. Both contributions are caused by the 
electron-electron interaction. However the kinetic contri- 
bution is present only out of equilibrium provided that 
the one-electron density of states is not constant near the 
Fermi energy. This contribution is proportional to the 
same combination of the one-electron energy distribution 



function as the inelastic relaxation rate and is thus inti- 
mately related with the relaxation. In the basis of exact 
many-body wavefunctions of the weakly interacting elec- 
tron gas the kinetic contribution is strictly off-diagonal 
in contrast to the thermodynamic contribution (persis- 
tent current) that depends only on the diagonal matrix 
elements of the current operator. The sign of the kinetic 
contribution is not fixed by the basic symmetry of the 
problem (orthogonal or symplectic) but depends on the 
nature of the energy dependence of the one-electron DoS. 

We are grateful to Igor Aleiner for critical comments 
and to B.L.Altshulcr for fruitful discussions clarifying the 
physical meaning of the kinetic contribution to the cur- 
rent. 
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